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ABSTRACT 



We describe deep (40 orbits) HST/STIS observations of the BALQSO PG 0946+301 and 
make them available to the community. These observations are the major part of a multi- 
wavelength campaign on this object aimed at determining the ionization equilibrium and 
abundances (lEA) in broad absorption line (BAL) QSOs. We present simple template fits to 
the entire data set, which yield firm identifications for more than two dozen BALs from 18 
ions and give lower limits for the ionic column densities. We find that the outflow's metalicity 
is consistent with being solar, while the abundance ratio of phosphorus to other metals is 
at least ten times solar. These findings are based on diagnostics that are not sensitive to 
saturation and partial covering effects in the BALs, which considerably weakened previous 
claims for enhanced metalicity. Ample evidence for these effects is seen in the spectrum. 
We also discuss several options for extracting tighter TEA constraints in future analyses, and 
present the significant temporal changes which are detected between these spectra and those 
taken by the HST/FOS in 1992. 

Subject headings: quasars: absorption lines — quasars: individual (PG 0946+301) 



1. INTRODUCTION 

Broad Absorption Line (BAL) QSOs are a spectacular manifestation of AGN out- 
flows. BALs are associated with prominent resonance lines such as Civ A1549, Si IV A1397, 
Nv A1240, and Lya A1215. They appear in about 10% of all quasars (Foltz et al. 1990) 
with typical velocity widths of ~ 10, 000 km s~^ (Weymann, Turnshek, & Christiansen 1985; 
Turnshek 1988) and terminal velocities of up to 50,000 km s^^. The small percentage of 
BALQSOs among quasars is generally interpreted as an orientation effect, and it is probable 
that the majority of quasars and other types of AGN harbor intrinsic outflows (Weymann 
et al. 1991). 

Establishing the ionization equilibrium and abundances (lEA) of the BAL material is a 
fundamental issue in quasar studies (Weymann, Turnshek, & Christiansen 1985; Wampler, 
Chugai & Petitjean 1995; Hamann 1996; Korista et al. 1996; Arav et al. 2001; de Kool et al. 
2001). Furthermore, determining the TEA is crucial for understanding the dynamics of the 
flows, especially radiative acceleration scenarios that are strongly coupled to the ionization 
equilibrium (Arav, Li & Begelman 1994; de Kool & Begelman 1995; Murray et al 1995; 
Arav 1996; Proga, Stone & Kallman 2000). Finally, the mass flux and kinetic luminosity 
associated with the flow cannot be constrained without a reliable lEA determination. 

Careful spectroscopic analysis has made it apparent that the BALs are often saturated 
while not black (Arav 1997; Telfer et al. 1998; Churchill et al. 1999; Arav et al. 1999b; de Kool 
et al. 2001). Support for this picture comes from spectropolarimetry studies (Cohen et al. 
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1995; Ogle et al. 1999; Brotherton et al. 2001). As a result, BAL ionic column densities (A^jon) 
determined in the traditional way (from equivalent width or direct conversion of the residual 
intensity into optical depth) can only serve as lower limits to the real column densities. 
Since BAL column densities are the foundation of any attempt to determine the TEA in the 
outflows, it became clear that more sophisticated analyses are essential for any progress in 
understanding the BAL phenomenon. To derive real Nion we have to account for saturation 
and partial covering factor in the BALs. In the UV, this approach is currently feasible only 
for exceptionally bright BALQSOs and still requires long integration times. After studying 
the available BALQSO data in the HST archive and hundreds of ground based spectra, we 
concluded that PG 0946+301 is the best candidate for such analysis. It has a substantially 
higher flux between 1200-2500 A (observed frame) than any other BALQSO that is suitable 
for lEA studies. Thus, we can obtain very high-quality data for a wide rest frame spectral 
region (500 - 1700 A), and the object suffers only minor contamination by Lya forest lines 
due to its low redshift {z = 1.223). Finally, the BALs of PG 0946+301 are typical for a 
high-ionization BALQSO. For example, the balnicity index of the CiV trough is 6300 km 
s~^, compared with 5000 km s~^ for the BALQSO sample of Weymann et al. (1991). 

It is for these reasons that many QSO researchers joined together to study PG 0946+301 
with deep, high-quality mult i- wavelength observations. The components of the campaign 
were: HST UV spectroscopy (described here), ASCA X-ray observations (Mathur et al. 
2000), FUSE UV spectroscopy (25 ksec data currently under analysis by the FUSE PI 
team), high- resolution optical spectroscopy and optical spectropolarimetry (some data al- 
ready taken, other observations planned). 

In this paper we describe the Space Telescope Imaging Spectrograph (STIS) observa- 
tions, present simple template fits for the BALs and discuss our preliminary lEA findings. 
We make the data (including some ground based spectra) available to the community in 
order to facilitate future analyses of this unique data set. In § 2 we elaborate on the data 
acquisition, reduction and wavelength calibration of the observations, and outline the tem- 
poral changes in the BALs detected by comparing the 1992 HST/FOS spectrum to the STIS 
data. We defer most of the temporal analysis to a future paper. In § 3 we describe our 
template fitting analysis for the entire spectrum, including the derivation of optical depth 
templates, line identifications and unidentified absorption features, and also discuss possible 
improvements over the template fitting technique. In § 4 we describe the lEA constraints 
obtained from the current analysis and discuss the over-abundance of phosphorus. In § 5 we 
summarize our results. 
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Fig. 1. — Composite STIS spectrum of PG 0946+301. Flux is measured in the observed 
frame. The bottom labels identify the ions responsible for the major absorption features. 
Some of the broad emission lines are marked above the spectrum. The dotted line designates 
zero intensity and the solid line above it is the noise spectrum. 
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2. OBSERVATIONS 

2.1. Data Acquisition and Reduction 

We observed PG 0946+301 with the Space Telescope Imaging Spectrograph (STIS) on 
the Hubble Space Telescope (HST) over a three-month period in early 2000. All of the spectra 
were obtained through a 52" x 0".2 slit to maximize throughput without a significant loss in 
spectral resolution. We used the MAMA detectors and G140L and G230L gratings to obtain 
full UV coverage at a resolution of 1.2 A and 3.2 A, respectively. We also used the CCD 
detectors and the G430M grating to obtain limited coverage in adjacent optical regions at a 
resolution of 0.56 A. The details of the observations are given in Table 1. 

Each G140L or G230L entry in Table 1 represents a five-orbit visit, since STIS MAMA 
observations are limited to five consecutive orbits due to constraints imposed by the South 
Atlantic Anomaly. To maximize the signal-to-noise of the spectra, we obtained exposures at 
five different locations separated by 0".3 along the slit (one location per orbit). We used this 
strategy to place the spectra at slightly different positions on the detector, to avoid possible 
fixed pattern noise not removed by the fiat-fielding of the spectra (Kaiser et al. 1998). For 
the G430M observations, we obtained two spectra per orbit (at the standard slit locations) 
to accumulate enough images to allow for accurate cosmic-ray rejection. 

We reduced the spectra using the IDL software developed at NASA's Goddard Space 
Flight Center for the STIS Instrument Definition Team (Lindler 1999). The individual 
spectra from each orbit were independently calibrated in wavelength and fiux, and resampled 
to a linear wavelength scale (retaining the same average dispersion). Careful examination of 
the spectra from each orbit reveals no discernible change in fiux over time or as a function 
of position on the detectors. We therefore averaged the spectra (weighted by exposure 
time) to obtain final versions of the G140L, G230L, and two G430M spectra. These were 
then combined into a single spectrum using Galactic absorption lines for a final wavelength 
calibration. The resultant spectrum and its noise are shown in Figure 1. 



2.2. Comparison of the EOS and STIS Data 

PG 0946+301 was observed by the HST Faint Object Spectrograph (FOS) in 1992. 
These data are still among the highest quality UV spectra of BALQSOs (complete analysis 
of these observations is found in Arav et al. 1999a, hereafter FOS99). Combined with the 
much higher S/N STIS data, we possess a unique data set for studying the temporal changes 
in the many BALs observed in the UV. In addition, we have five epochs of the C iv trough 
obtained from ground based telescopes. We defer a full analysis of the observed temporal 
changes to a separate paper. Here we briefiy describe the main changes and concentrate on 
the ones that affect our model fits most. 
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Fig. 2. — The Civ spectral region in PG 0946+301 at two different epochs. Histogram: 
STIS data from March 2000; dotted hne: FOS and Lick data from 1992. The two settings 
of the G430M STIS grating are shown (rebinned by 10 pixels), these overlap between 3165 - 
3300 A. In the main trough (3260 - 3400 A) there is a significant deepening of the absorption 
in the 2000 data and a new high velocity trough is seen between 3130 - 3250 A. We note 
that the differences are very significant compared to the noise level and that the STIS data 
are in remarkable agreement with Lick data from the same epoch (not shown here). 
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In Figure 2 we show the Civ spectral region at two epochs separated by 8 years (~3.5 
years in the rest frame of the object). The 1992 spectrum is a combination of the HST 
data from the FOS G270H grating (A < 3250 A, with a four pixel rebinning), and Lick 
data (A > 3250 A). The 2000 observations consist of two setting of the STIS G430L grating 
rebinned by 10 pixels (see Table 1 for exact dates and exposure times). Our Lick 2000 data 
(not shown here) are in remarkable agreement with the STIS G430L data. Inspection of 
Figure 2 reveals several significant changes between the 1992 and 2000 spectra. 

1. Significant deepening of the main trough has occured between 3260 - 3400 A (observed 
frame). For comparison the S/N level at the deepest part of the trough is 6.5 for the 
Lick 1992 data and 7 for the binned STIS data. 

2. A substantial decrease in flux is seen between 3280 - 3300 A (~ 25%). A similar change 
is also easily identified in the N v and O vi BALs. 

3. A new shallow high velocity trough is seen in the STIS data between 3130 - 3250 
A (-17,200 to -28,400 km s^^), with a maximum depth of 20%. Several other high 
ionization BALs show this feature (most evident in N v and O Vl) . The high velocity 
trough is not present in low ionization species (Si IV, Cm, Nlll...). 

4. The Si IV absorption depth in the broad low velocity trough (-4000 to -7000 km s~^) 
has increased significantly (not shown in the figure), and this can also be seen in other 
low-ionization species (e.g. Olll, Nlll). 

Both independent settings of the G430L grating show the new high velocity BAL as well 
as the decrease in flux that is seen between 3280 - 3300A. We also note that while the Si iv 
BEL is 20% stronger in the STIS data, the Civ BEL has not changed appreciably (based 
on the Lick 2000 data). 



3. TEMPLATE FITTING ANALYSIS 

The simplest way of analyzing an absorption line spectrum containing many broad 
overlapping lines with a complex velocity structure is the template fitting approach (Korista 
et al. 1992; FOS99). This method uses the apparent optical depth (defined as r = —ln{Ir), 
where Ir is the residual intensity seen in the trough, for more details see FOS99) profile 
of a given ionic transition as a template to fit all other transitions, where only a simple 
multiplicative scaling is allowed. The fitting process finds the optimum scaling values which 
minimize the difference between the data and a model spectrum based on using scaled 
versions of the template for all plausible transitions. We did this in two ways. First, we 
used standard x^ statistics, i.e., minimized the square of the differences between the model 
and the data, weighted by the square of the Poisson error. However, because the S/N ratio 
is very high everywhere, the deviations between the model and the observations are almost 
entirely due to uncertainties regarding the effective continuum (defined below) and covering 



factor effects. Therefore, we also fit the spectrum using a constant error rather than the 
Poisson error at each wavelength point. In practice using a constant error or the statistical 
error changes the parameters derived from the fit by less than 0.1 dex so that the exact 
minimization method used is irrelevant. 

The fitted scaling factors are translated to ionic column densities (A^jon) using standard 
techniques (see below). Two assumptions must hold in order that the fitted N.ion be good 
approximations for the real ones. First, the absorbing material covers the emission source 
completely and uniformly, and scattered-photons do not contribute appreciably to the resid- 
ual intensity in the trough (see discussions in Korista et al. 1992; Arav 1997). Second, the 
ionization equilibrium does not vary as a function of velocity. As we will show, the first as- 
sumption is not a good approximation for the BALs in PG 0946-1-301. In this case we must 
treat the fitted Nion as lower limits for the real ones (FOS99, Arav et al. 1999b, de Kool et 
al. 2001). 



3.1. Methodology 

Using the assumptions described above, the optical depth at a given wavelength A due 
to a transition between levels / and u of species n can be expressed as 



Tnlu{^) — A/u// 



Tre^ ^ ^ Xni dNo 



lu' 



rrieC Xq dv 






Here A/^ and //„ are the rest wavelength and oscillator strength of the transition, ^ the 
density ratio between ions of species n in energy state / and some reference species for 
which the distribution of column density over velocity ^^ is known, and v{\iu -^ A) is the 
velocity that shifts a line with rest wavelength A/^ to the wavelength A. The distribution of 
column density over velocity ^^ is usually derived from the optical depth in a transition 
corresponding to a line that is not blended with other lines: 

dNo rrieC 1 



dv 7re2 Ao/o 



rM (2) 



It can be derived directly if the transition is a singlet (e.g. de Kool et al. 2001), or by a 
deconvolution if the transition is a multiplet (e.g., FOS99). The line profile To(f) of the 
reference transition thus serves as a template for all other absorption lines. 

Once the template is known, a model of the total optical depth as a function of wave- 
length is constructed based on a first guess of the ratio of the column density of each species 
relative to that of the reference species, ^^: 

^(^) = EEE^"'«(^) (3) 

n I u>l 
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Using the derived r(A) and a model for the unabsorbed continuum, a model spectrum is 
computed that can be compared with the observed spectrum. Using an iterative procedure, 
the values of ^ are then adjusted until the x^ is minimized. The value of ^ of the best fit 
yields the ratio of the column density of ion n in energy state / to the column density of the 
reference species. In the current analysis we assume that the ground-term level populations 
of the ions we consider are in LTE, at a temperature much higher than the level energy, so 
that the relative population of each energy level in the ground-term is proportional to its 
statistical weight. 



3.2. The Velocity Templates 

As in our previous study of PG 0946-1-301 we used templates derived from two reference 
transitions. Civ A1548.20 is used for lines from high ionization species, and Si IV A1393.76 
for lines from low ionization species, both lines are the blue component of their respective 
doublets. In Figure 3 we present the C IV and Si IV optical depth as a function of velocity 
derived from the STIS spectrum of PG 0946-1-301. The template for Civ was obtained with 
a regularization method similar to that described in FOS99. Because of the much higher 
resolution in the observed Civ line, the regularization used in FOS99 which suppressed 
fluctuations in neighboring points did not give satisfactory results. Instead, we used an 
alternative form of the regularization matrix H, in the notation of FOS99: hi^i = 2, /ij,j+s = 
—1, hi^i-s = — 1 and all other elements of i/ equal to zero. This form suppresses oscillations on 
the scale of the doublet separation but does not smooth the derived profile over neighboring 
points. The value of the regularization parameter used was 0.1. We derived the Si iv template 
without regularization, since oscillations on the scale of the doublet are not apparent in the 
plain deconvolution result. 

The transition between data from the G230L grating and the G430M grating, and the 
accompanying change in spectral resolution, occurs in the middle of the Si IV line. This 
accounts for the smooth appearance of the high velocity part of the Si IV template. There is 
no evidence of the new high velocity trough in the Si iv absorption line. The Si iv template 
was corrected for the presence of intervening C IV absorption at 1330 A in the QSO restframe 
by removing this line from the spectrum before the deconvolution. We note that some of 
the absorption in the range -11,000 to -13,000 km s~^ may be due to CiV from another 
intervening system that can be identified from its Lya absorption (component 9 in Table 
2 in FOS99). There is some evidence for this being the case from the fact that the Si IV 
template over-predicts absorption in other low ionization species (eg, O ill, S iv) at these 
velocities. 
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Table 1. HST/STIS Observations of PG0946+031 



Grating 


Range 


Exposure 


Date 




(A) 


(sec) 


(UT) 


G140L 


1140- 1715 


12,960 


2000 February 26 


G140L 


1140- 1715 


12,960 


2000 March 28 


G140L 


1140- 1715 


12,960 


2000 March 30 


G140L 


1140- 1715 


12,960 


2000 April 4 


G140L 


1140- 1715 


12,960 


2000 May 11 


G430M 


3022 - 3306 


5,274 


2000 April 12 


G430M 


3163 - 3447 


7,250 


2000 April 12 


G230L 


1640 - 3148 


12,960 


2000 April 23 


G230L 


1640 - 3148 


12,960 


2000 May 2 




-3x10 



■2x10 -1x10 

velocity [km/sec] 



Fig. 3. — The apparent optical depth as a function of velocity for the blue components of 
the Civ A1549, (thin line) and Si IV A1397 (thick line) resonance doublets. 
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3.3. The Effective Continuum 

In order to make measurements or fit the observed BALs we must use a model for 
the unabsorbed emission (continuum plus BELs, hereafter "the effective continuum" ) of the 
quasar. At wavelengths longer than ~ 1000 A, the spectrum of PG 0946+301 is very close to 
the non-BALQSO composite of Weymann et al. (1991), which we therefore used as our long 
wavelength effective continuum with only small modifications to some emission line strengths 
(as in FOS99). Shortwards of 830 A we found only two narrow regions (624 - 631 A and 
567 - 571 A) free from BAL features. This complicates the task of determining the effective 
continuum based on the data at hand. Unfortunately, very little is known observationally 
about intrinsic QSO spectra in the range 500 - 900 A. The QSO composite spectrum of Zheng 
et al. (1998) does cover these wavelengths, but very few QSOs contribute to the composite 
in this wavelength range. Furthermore, the signal to noise ratio is so low that the use of this 
composite as a continuum would introduce too much noise to the analysis. 

With these uncertainties in mind, we used a similar approach to the one we employed in 
the analysis of the HST FOS data. Three effective continua are employed (shown in Fig. 4). 
The first is a fixed continuum which is very similar to the fixed continuum in FOS99 (for 
details see § 5.2 of FOS99). We lowered the fiat continuum shortwards of 900 A by 5% and 
added two BEL features, one for O v A630 and one for Ne v A571. Both BELs are required by 
features evident in the STIS data (the quality of the FOS G160L data was too low to justify 
the addition of these BELs). Emission features of roughly the same shape and wavelengths 
are also seen in a smoothed version of the quasar composite spectrum produced by Zheng 
et al. (1998). 

The second and third effective continua are different from the fixed continuum for wave- 
lengths below 1060 A. In this range, the underlying continuum fiux Fc was assumed to have 
the functional form 

F, = Fo + Fi(^)W (4) 

where Fq, /5, Xcut and Fi are free parameters with the constraint that the continuum matches 
the long wavelength one at 1060 A. The values of the continuum parameters were determined 
by optimizing them simultaneously with the column densities of all species (^f^) during 
the template fitting process. Since the continuum parameters were allowed to vary in the 
fitting process, we refer to the resulting continuum as a fioating continuum. The difference 
between the two floating continua is that one is derived from a template flt that only includes 
resonance lines as opacity sources, whereas the other also allows for a Lyman edge opacity 
(as a free parameter). The equivalent widths of the model BELs were flxed during the fltting 
process. Finally, most of the known intervening absorption line systems were modeled with 
single Gaussians with flxed equivalent widths to minimize possible confusion with BAL 
structures. 
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3.4. Results of the Template Fits 

The column densities resulting from the template fits are summarized in Table 3, and 
the fit itself (for the fixed continuum) is illustrated in Figure 5. Before discussing these 
results, we re-emphasize that column densities obtained from template fitting are in reality 
only lower limits because of partial covering effects. For some species we will also quote 
an upper limit on the column density. Such an upper limit can be derived if the species in 
question has at least one line with a much lower oscillator strength than the others, which 
is not detected. These upper limits are necessarily somewhat subjective since they depend 
on the shape of the effective continuum and possible blending with other lines. By stating 
that a line is not present we are in fact making a judgment as to what variations in the 
effective continuum would be reasonable. Therefore, the upper limits are derived in practice 
by manually changing the column density in the fit and comparing the model spectrum and 
the observed spectrum at the position of the weak line to see what column density could be 
present without conflicting with the observations. 

When inspecting Table 3, the first thing to note is that the derived column densities are 
insensitive to the exact choice of effective continuum. Even changes of 30% in the continuum 
level affect the column densities of the stronger lines by a similar amount (i.e., 0.1 dex in the 
table), while the effects of non-black saturation can increase the inferred column densities 
by a factor of three or more. 

Here follows a summary of the individual results: 

Hydrogen. The Lya BAL is not very deep, leading to a low value for the column density 
(\og{NHi) = 15.4) derived from it. If the fit includes a Lyman edge due to the BAL outflow, 
the best flt Lyman edge corresponds to a column density of log(A^Hi) = 16.4. The neutral 
hydrogen column density from the Lyman edge flt is still consistent with the upper limit 
derived from the absence of a clear signature for Ljf3 and Ly7 BAL, log(A'^^i) < 16.7. 
In our ionization discussion (§ 4), we use a somewhat more conservative upper limit of 
log{NHi) < 17 to allow for partial covering effects. 

Helium. The upper limit quoted for He I is derived from the absence of a clear signature 
of the 584 A resonance line. This is illustrated in Figure 6a. 

Carbon. We do not detect any absorption feature associated with a Cll A1335 BAL. 
The upper limit for this line provides important ionization equilibrium constraints (see § 4). 
A C III BAL is clearly detected and its deepest component is not severely blended with other 
BALs, so that the apparent column density is well deflned. The maximum optical depth is 
very similar to that of other low ionization species (Hi, Nlll, Siiv), which strongly suggests 
that the low ionization BALs are essentially saturated and derive their shapes primarily 
from partial covering. The C IV BAL is one of the templates, so no new information can be 
derived from its template flts. 
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Fig. 5. — Expanded view of a fit and the fixed effective continuum used in producing it. Most 
BAL features are identified by the absorbing ion. The short vertical fines mark intervening 
and Galactic absorption lines, which we superimpose on the effective continuum to account 
for non-BAL absorption features. 
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Nitrogen. N ill gives rise to several lines in different parts of the spectrum. The value 
for the NlII column density is driven mainly by the strongest lines around 685 A, the other 
ones (991 A and 764 A) being both weaker and blended with stronger lines from other 
species. The column density for Niv is not very well determined since its BAL is a part 
of the large blended trough near 750 A. The Nv BAL is very deep, even though it falls 
under the Lya emission line. This contrasts with the saturation effects in the low ionization 
species mentioned above, which set in at a much higher flux level and imply a low covering 
factor. The higher ionization BALs appear to cover a significantly larger part of the source 
than the low ionization BALs, including most of the Lya BEL region.. This suggests a 
physical picture where at a given velocity interval there is lower density extended material, 
surrounding denser lower ionization material with smaller geometrical extension, (see FOS99 
for a full discussion). The new shallow high velocity component is deeper in Nv than 
expected based on scaling the CiV template, as illustrated by the fact that the flux in the 
fitted model significantly exceeds the observed flux at the wavelengths where this absorption 
occurs. Note that this is not likely to be a result of the uncertainty in the continuum level 
since a comparison between the current STIS observation and the 1992 FOS observation 
clearly shows the emergence of the new absorption component. The fact that the high 
velocity component is deeper in Nv than in CiV suggests that it is highly ionized, which is 
supported by the strength of this component in O VI and S VI. 

Oxygen. Oxygen is a unique element in this analysis since we have relatively unblended 
BALs for four consecutive ionization stages. O ill has two multiplets at 703 A and 835 A. 
The O III BALs are reasonably well fit by the model. The maximum optical depth in O ill 
is higher than in the other low ionization species, which is consistent with the trend that 
ions with higher ionization potential (IP) tend to have a larger covering factor (The O ill IP 
is almost midway between the Si IV and CiV IPs). Unfortunately, there are no Olll lines 
that allow us to determine an upper limit for the column density. O iv has three groups of 
lines, and the fit is mainly driven by the one with the lowest oscillator strength around 609 
A. The other two BALs (554 A and 789 A) are part of the large troughs around 550 A and 
750 A. If the 609 A BAL is well fit, the other two BALS are predicted to be deeper than 
observed, giving rise to the very large Tmax- Because this occurs close to the bottom of the 
trough, the difference between model and fit does not contribute as much to the overall x^- 
Again, this demonstrates the effects of partial covering. The O V BAL is interesting in that it 
clearly demonstrates how complex the real BAL formation is compared to our simple model. 
Ov A630 is a singlet and the main part of the absorption (from -4000 km s~^ to -11,000 
km s~^) is not blended with any other BALs. If our assumptions were correct, we would 
expect the shape of the O v BAL to correspond well to the C iv velocity template. In reality, 
however the shape is quite different, with the main velocity component in the template at 
-10,000 km s~^ being less deep than the broad component between -4000 km s"^ and -9000 
km s~^. A possible explanation for this behavior is that in this part of the spectrum the 
source size is varying rapidly as a function of wavelength, perhaps because of the presence 
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of partially uncovered emission lines. We consider it more likely that the problem lies with 
continuum/source size problems around 600 A than with the template, since many other 
high ionization lines, including O VI, were fit very well with the C iv template. O vi also 
clearly shows the new shallow high velocity component. 

Neon. Absorption lines from several ionization stages of neon are expected to occur in 
the spectrum of PG 0946+301, and the template fitting procedure finds very large column 
densities for all of them. However, the neon lines occur only within the two large blended 
troughs around 550 A and 750 A so that their apparent column densities are less reliable 
than those of most other ions. The exception is Ne v, its A571 multiplet is the main source 
of opacity for the red edge of the 555 A trough, and therefore its fitted value is quite reliable. 
The combination of the uncertainty in the real optical depth in the neon lines with their 
low oscillator strengths leads to large values for the column density that are not necessarily 
realistic. However, establishing the existence (or lack) of a Nevill BAL is of particular 
importance in the context of ionization equilibrium constraints (see § 4). We therefore 
conducted a tailored fitting experiment to assess its existence. We fit only the 700 - 845 A 
region and ignored the correlation with other lines from the same ion found in different parts 
of the spectrum. In this way we circumvent some of the non-physical restrictions associated 
with ignoring saturation and partial covering inherent in our fitting process, and obtain the 
best local fit possible. Figure 6b shows the data, a fit including Nevill and a fit excluding 
it. The large improvement caused by the inclusion of Ne vill is clearly seen in the 768 - 777 
A 753 - 757 A and 730 - 740 A regions. We conclude that the presence of an appreciable 
Nevill BAL (log(A^jvevin) ^ 16.5) in PG 0946+301 is highly probable. 

Magnesium. In our analysis of the FOS data we discussed the probable identification 
of the Mgx AA 609,625 BAL, which was based on fits to low-quality data obtained by the 
G160L grating. Our vastly superior quality STIS G140L data do not show any features that 
can be attributed to this BAL. The excellent fit for the 580 - 605 A region shown in Figure 
5 does not include any Mg X contribution. However, the strong presence of BALs from O IV 
and O V do not allow the determination of an interesting upper limit for this ion, the value 
we obtain for the upper limit is 1-3 xlO^^ cm~^. 

Silicon. We do not detect any absorption feature associated with a Si ll A1263 BAL. The 
upper limit for this line supplies similar ionization equilibrium constraints to those available 
from Cll (see § 4). A possible Sim A1206 BAL would be blended with the Lya and Nv 
BALs, so that very tight constraints are not possible. There is no evidence of the main 
component at -10,000 km s~^ in the spectrum, leading to the upper limit in Table 3. Like 
C IV, Si IV is one of our template species, therefore no new information is obtained from the 
fit. 

Phosphorus. As in previous investigations, we find good evidence for the presence of 
a Pv BAL. In spite of its high ionization potential, the profile of the Pv BAL fits much 
better with the Si IV template than with the C iv template. This would be consistent with 
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a picture in which the BAL region contains a component with a very high column density 
and a small covering factor, which is responsible both for lines of lower ionization species 
and for lines from very low abundance species (see FOS99 for a full discussion). A possible 
P IV A950 BAL would occur in a highly blended region of the spectrum, so that only a 
relatively uninteresting upper limit can be derived. 

Sulphur. S III has a large number of lines that are predicted to form very broad BALs 
without sharp features. Because of this, identification of S ill features is not very secure, 
and S III is only marginally detected in the spectrum, as is illustrated in Figure 6c. S IV is 
very clearly present, and the STIS spectrum resolves the uncertainties associated with the 
S IV hues in the FOS spectrum of PG 0946+301 (see FOS99). The depths of the strongest 
S IV lines is similar to that of the low-ionization template. S iv gives rise to several lines 
with a range of oscillator strengths. This allows a determination of an upper limit on the 
Siv column density of logA^(S'lv) < 15.8 (based on Siv AA809,815), which makes Siv the 
species with the most tightly constrained column density in our analysis. Sv on the other 
hand is very poorly constrained, and the global template fits yield a negligible column density 
for it. The reason for this is that the S v line almost coincides with one of the stronger OlV 
lines at 790 A. As mentioned above, the O IV fit is driven by a very weak line that already 
over-predicts O iv absorption at the position of the S v BAL, so that the S V column density 
is driven to zero. Since S v has only one line, this can occur without incurring penalties for a 
bad fit in another part of the spectrum. The value quoted in Table 2 is derived from a local 
fit (700 - 850 A), which suffers less from the effects described above. S VI is the strongest and 
deepest sulphur BAL, in agreement with the hypothesis that higher ionization species have 
a larger covering factor, although the BAL definitely has a profile that is slightly different 
from that of the other high ionization BALs. The shallow high velocity component is evident 
in S VI, where it is deeper than that in C IV. 

Argon. Resonance lines of Ariv, Ar v, Ar vi and Arvil occur in the wavelength range 
covered by the STIS spectrum of PG 0946+301. None of these lines are detected. For 
Ariv, Arv and Arvi this leads to upper limits to their column density ~ 2 x 10^^ cm~^ 
because these species have lines with fairly small oscillator strengths. These limits do not 
provide interesting constraints. The situation is different for Ar vil which has a strong singlet 
(oscillator strength of 1.24) at 586 A. Figure 6d shows a fit which includes the Arvil line at 
the level of 1.3 x 10^^ cm^^. However, the fit is unphysical since at the expected position 
of maximum Ar VII optical depth (567 A) we find a much higher residual intensity than the 
fitted value. (We note that the adjacent narrow absorption feature seen in the data just 
longward of this position is due to Galactic Sill A1260). In all the well identified BALs the 
deepest predicted feature is fitted much better. In addition the Arvil fit does a very poor 
job in reproducing most of the surrounding spectral features. Based on these arguments, 
10^^ cm~^ is a conservative upper limit for the Arvil column density. 

Iron. The spectrum was searched for the presence of BALs from Fe ill and Fe iv. Fe ill 
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Fig. 6. — Fit tests used to determine the existence of marginal, or heavily blended species 
(see text for discussion). The observed spectrum is shown as histograms, the thick and thin 
lines represent models with and without the ion considered, respectively. 
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has approximately 130 transitions in the wavelength range considered, and if present, would 
give rise to relatively broad blended BALs, similar to the case of S ill. The strongest lines 
occur around 860 A, and the absence of any detectable features corresponding to these leads 
to logA^(Felll) < 15.3. Feiv has a resonance multiplet centered around 525 A. Including 
FelV does improve the fit at the very shortest wavelengths, but there is no clear feature 
corresponding to it. 

Unidentified absorption features. We found several absorption features in the STIS data 
for which we do not have good identification with either expected BALs, or with intervening 
and Galactic absorption lines. Most prominent are the features at (observed wavelength): 
1526 A (r ~1), 1484 A (r -0.3), 1463 A (r -0.3), 1277 A (r ~1). We note that the feature 
at 1526 A is too deep and broad to be fully explained by Galactic Sill A1526 line, since it 
should not be stronger than the observed Sill A1260 feature. 
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Table 2. Fitted Column densities 



Ion 



Lines 



Template '^ Nion 



Tmax ^[Tmax) ^^ ion ' ^^ ion 



Hi 


1215.67 


0.416 


A 


15.4 


0.82 


1215.67 


- 


- 


Ly lim 


911.76 


- 


A 


- 


- 


- 


- 


16.4 


He I 


584.33 


0.285 


A 


<15.5 


- 


- 


- 


- 


Cll 


1335 


0.127 


B 


<14.7 


- 


- 


- 


- 


Cm 


977.02 


0.767 


B 


15.1 


0.82 


977.02 


15.2 


15.2 


Civ 


1549 


0.285 


A 


16.1 


3.34 


1549 


- 


- 


Nm 


991 

764 
686 


0.122 
0.082 
0.402 


B 


15.6 


0.84 


686 


15.6 


15.6 


Niv 


764.15 


0.616 


A 


16.1: 


4.22 


764.15 


16.1: 


16.0; 


Nv 


1240 


0.235 


A 


16.2 


2.88 


1240 


- 


- 


Qui 


834 
703 


0.107 
0.137 


B 


16.3 


1.31 


834 


16.3 


16.2 


Oiv 


789 
609 
554 


0.110 
0.067 
0.335 


A 


16.9 


8.83 


554 


16.8 


16.8 


Ov 


629.73 


0.515 


A 


16.2 


3.30 


629.73 


16.2 


16.2 


Ovi 


1034 


0.199 


A 


16.6 


5.04 


1034 


16.6 


16.6 


Neiv 


543 


0.233 


A 


16.6: 


2.68 


543 


16.6: 


16.6: 


Nev 


571 


0.092 


A 


16.8 


1.68 


571 


16.8 


16.6 


Nevi 


561 


0.090 


A 


16.8:: 


1.56 


561 


16.8:: 


16.8: 


Neviii 


774 


0.302 


A 


16.5: 


2.36 


774 


16.5: 


16.5: 


Sin 


1263 


1.18 


B 


<14.2 


- 


- 


- 


- 


Sim 


1206.50 


1.68 


B 


<14.3 


- 


- 


- 


- 


Si IV 


1397 


0.784 


B 


15.1 


0.86 


1397 


- 


- 


Piv 


950.66 


1.16 


B 


<15.0 


- 


- 


- 


- 


Pv 


1121 


0.701 


B 


15.0 


0.36 


1121 


- 


- 


Siii 


680 


1.64 


B 


14.9:: 


0.41 


680 


14.6:: 


14.6: 


Siv 


814 
750 
660 


0.103 
0.749 
1.17 


B 


15.3 


1.06 


660 


15.2 


15.1 


Sv 


786.47 


1.42 


A 


15.2™ 


0.70 


786.47 


- 


- 


Svi 


937 


0.665 


A 


15.8 


1.81 


937 


15.8 


15.8 


Arvii 


585.75 


1.24 


A 


<15.1 


- 


- 


- 


- 


Feiii 


861 


0.131 


B 


<15.3 


- 


- 


- 


- 



^Lines that play a part in our analysis. Multiplets are rounded to the nearest A. 

Total oscillator strength of the niultiplet. 
'^A refers to the C IV template, B to the Si IV template. 
'^\og(Nic,n) in units of cm^'^. 
°Fixed continuum. 

The maximum optical depth due to the ion. 

^Wavelength (A) of the multiplet for which Tmax occurs. 

'^Floating continuum (lines below 1060 A only). 

''Floating continuum including Ly edge (lines below 1060 A only). 

™The S V BAL is strongly blended with Nevm, Niv and Oiv BALs. Therefore, we only used a regional fit (700-845 A) with 
the fixed continuum to constrain its column density. 

(:) spectral features from this species definitely present, but column density uncertain. (::) spectral features from this species 
only marginally detected. 
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3.5. Beyond Template Fitting: Deriving and Constraining Real Nion 

As discussed in the introduction, a main goal for further analysis of the PG 0946+301 
STIS data is to derive real Nion for the BALs or at least useful constraints beyond the simple 
template lower limits. While not attempting such a full analysis here, we outline several 
promising approaches as well as potential difficulties for doing so. 

An obvious place to start is to examine the departure of the template fits from the data. 
In cases where we have BALs from two or more lines from the same ion, these departures 
often indicate non-black saturation and may be used to determine the covering factor of the 
absorber. Two good examples are the O IV lines A609 and A554 with an oscillator strength 
ratio of 1:5. Inspection of our fits clearly shows that the A554 line is expected to be black 
if the A609 is well fit. Since the A554 trough is not black, we conclude that this is a case 
of non-black saturation where the covering factor can probably be measured directly from 
the shape of the A554 trough. A similar case is found for Nlll A991 and A686. In order to 
estimate the degree of saturation, a comparison can be made between global fits for the N ill 
A991 and A686 lines and fits that assume these lines are independent. As a refinement one 
can try fits where individual multiplet components are saturated. Probably the easiest Nion 
to measure is that of Siv. We detect four Siv multiplets at 1070 A, 814 A, 750 A and 660 
A, with oscillator strength ratio of 1:2:15:24. Modeling of the troughs associated with these 
multiplets should yield an even more accurate Nion for S iv than we have currently. 

A large volume of work (e.g., Arav 1997; Barlow 1997; FOS99; Churchill et al. 1999; 
Arav et al. 1999b; de Kool et al. 2001) shows that both the level of saturation and the 
covering factor can be strongly dependent upon the velocity. Therefore, the extraction of 
Nion and the following determination of the lEA in the fiow should be done for as many 
velocity components as possible. Ideally, it would be best to work with the column density 
as a function of velocity instead of the velocity integrated column density. 

The task of extracting Nion is complicated by the following difficulties: 1) Most doublet 
and multiplet components are blended together since the velocity width of the BAL is much 
larger than the separation between the different components of the same line. If the source is 
completely covered and there is no scattered light component at the bottom of the troughs, 
an optical depth solution can be extracted from blended components (Junkkarinen, Burbidge 
& Smith 1983, FOS99). However, this is obviously not the case for most BALs present in the 
spectrum of PG 0946+301. 2) The considerable width of the BALs causes blending between 
different lines. The most extreme example is the deep trough between 720 - 770 A, which is 
a blend of at least four strong BALs and probably a few additional weaker ones. 3) The new 
high-velocity trough associated with the high ionization species (Civ, Nv, Ovi...) worsens 
the BAL-blending problem. Although not very deep, the high-velocity trough extends over 
a velocity width comparable with that of the main deep trough. Therefore, the combined 
velocity width of the BAL system in a given high ionization line has doubled compared with 
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the situation in 1999, which worsens the BAL blending. 



4. CONSTRAINTS ON THE lEA IN PG 0946+301 

The template fitting results yield important constraints for the ionization equilibrium 
and abundances (lEA) in PG 0946+301. All the available lower limits on Ajon, as well 
as upper limits based on non detections, can be used to constrain the lEA. The following 
findings are particularly useful in doing so: 

1. There is no clear signature for either Ly f3 and Ly7 BALs, or a BAL Lyman-limit, which 
yield the upper limit Nhi ^ 10^^ cm~^. Our template fit for the Lya BAL gives the 
lower limit of Nhi k, S x 10^^ cm"^. 

2. Lower limits for the ionization parameter can be derived from the absence of neutral 
helium (inferred from the non detection of a He I A584 BAL) and singly ionized metals 
(Cii, Nil, Sin...) 

3. A departure from solar abundances ratio is implied by the presence of P v combined 
with the absence of Arvil, the Hi upper limit and the Siv constraints. 

In this paper we concentrate on constraining simple-slab photoionization models (using 
the photoionization code CLOUDY; Ferland 1996). Such models are commonly used in 
the study of quasar outfiows (e.g., Weymann, Turnshek, & Christiansen 1985; Arav & Li 
1994 Hamann 1996; Crenshaw & Kraemer 1999; de Kool et al. 2001) and assume that 
the absorber consists of a constant density slab irradiated by an ionizing continuum. The 
two main parameters in these models are the thickness of the slab as measured by the 
total hydrogen column density {Nu) and the ionization parameter U (defined as the ratio 
of number densities between hydrogen ionizing photons and hydrogen in all forms). The 
spectral shape of the ionizing continuum is also important. For this analysis we used two 
specific spectra: The first is the standard Mathews-Ferland (MF) AGN spectrum (Mathews 
& Ferland 1987). The MF spectrum is very hard between 1-4 Ry, on the average f,y oc z/""'^ 
over that range (the so called big blue bump). In PG 0946+301 we do not see such spectral 
characteristics. Instead, we observe f,y oc z/~^ between 1-1.5 Ry with a hint of an even softer 
slope thereafter (see Fig. 4). We therefore constructed a more realistic ionizing-continuum 
model for PG 0946+301. Below 1 Ry and above 40 Ry we match the spectral shape (and 
relative intensities) of the MF spectrum, between 1-40 Ry we use an /^ oc z/~^ slope. This 
shape is consistent with the observed spectrum while keeping the shape of the far UV - soft 
X-ray spectral region as simple as possible. We note that the expected Galactic extinction 
in the direction of PG 0946+301 is negligible (E[B-V] = 0.018 mags., Schlegel, Finkbeiner & 
Davis 1998), and therefore should not cause an appreciable correction for the intrinsic UV 
spectral slope. The results described below were derived using this latter spectral shape. For 



-23- 



comparison we discuss the quantitative differences which arise from using the MF spectrum 
at the end of this section. 

A particularly useful way to constrain such models is to plot curves of constant Ni^n 
on the plane of Nh vs. U . On such a plot any type of information regarding the ionic 
column densities can be used to constrain the allowed parameter space. Lower limits (derived 
from fits and/or apparent optical depth measurements) exclude the area below the constant 
Nion curve, Ni^n upper-limits (from non-detections) exclude the area above their associated 
constant-iVjon curve. In Figure 7 we show the most significant Nion constraints extracted 
from the STIS data where we use the fi, oc z/~^ spectrum described above as the incident 
ionizing continuum and assume solar abundances. For Hi we have both a lower limit from 
our fit to the Lya BAL and an upper limit from the absence of clear signature for either 
Ly/9 and Ly7 BALs, or a BAL Lyman-limit. In our Nh vs. U presentation, the combined 
constraints restrict the allowed parameter space to the diagonal line-shaded band. A similar 
situation arises for Siv where we derive a lower limit on A^^iv from the strong multiplet at 
660 A and an upper limit from the weaker multiplets at 1070 A and 814 A. The combination 
of both constraints gives the dot-shaded area. 

Several important conclusions can be drawn from Figure 7. The combined constraints 
from Hi, Siv, Cii, Arvil as well as from other available CNO ions (which for clarity's sake 
are not shown here), suggest that the absorber is characterized by 20.3 < log(A''^) < 20.8 
and —1.7 < log(f/) < —1.2. This narrow range in the Ng/U plane is fully consistent with 
23 out of our 27 ionic column-density constrains, without assuming departure from solar 
abundances. However, the limits on the allowed metalicity range are rather wide, reflecting 
the uncertainty in the Hi BAL column density (see § 3.4). We find that a simple scaling 
of the metalicity in our models allows for solutions consistent with the H I constraints for 
0.7Zq < ^0946 < IOZq. The total column density of the absorber is inversely proportional 
to its metalicity value. 

There are two constraints which are inconsistent with this narrow range in the Nu/U 
plane. The firm detection of P V (and even more so in the FOS data, see Junkkarinen et al. 
1997) yields the associated lower limit on Figure 7, which is incompatible with the allowed 
parameter space we derived above. The simplest way to resolve this apparent contradiction 
is to assume that phosphorus is over-abundant by a factor of 10-30 compared to its solar 
abundance. All else equal, this will lower the P v curve to the point where it intersects with 
the allowed parameter space. Further discussion of the inferred phosphorus over- abundance 
is found in § 4.1 and in the Appendix. 

Another inconsistency is evident from the Ne Vlll curve, for which the allowed parameter 
space is strongly excluded by most other metal constraints. Here there are two possibilities. 
If the Ne VIII measurement is robust, our naive assumption of a single U value (or a narrow 
range of U) must be wrong. A second option is to question the Ne vill detection since the 
Ne VIII BAL resides in the middle of the large trough and is blended with strong O IV, S V 
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and N IV BALs, as well as from a few weaker ones. However, as we show in § 3.4 and in Figure 
6 it is very probable that the wide absorption trough around 750 A contains a considerable 
Ne VIII BAL contribution. The constraints from S VI and Ne V are also marginally above 
the allowed area in the N„/U plane (by ~ 0.2 dex). Since both ions are of relatively high 
ionization potential, these results support the hypothesis of a larger range of higher U values 
in the outflow, which is needed to explain the inferred NeVlll value. We also point out 
that the uncertainties in the continuum's shape increase at higher energies (since our last 
data point lies near 1.8 Ry). Therefore, in the vicinity of these high ionization potentials 
the extrapolated flux level is rather uncertain, which may contribute to the inconsistency 
between the model results and the observed constraints for these ions. 

How relevant are these constraints to the actual PG 0946+301 outflow? Our analysis 
of the FOS data (FOS99) suggests a range in the ionization parameter at a given velocity. 
Furthermore, it shows that the velocity dependence of the apparent optical depth is different 
for low ionization BALs than for high ionization ones (some of this is accounted for by 
using different templates for high vs. low ionization species). We argue that these effects 
can only modestly change our overall lEA conclusions, but there might be some minor 
outflow components with differing constraints which are missed by ignoring the velocity 
dependencies. Several arguments support this assertion. The P V relative ionic fraction is at 
a maximum inside the U range already derived from the other metals/hydrogen constraints, 
— 1.7 < log(f/) < —1.2. Therefore, a wide range in U will only increase the inferred over- 
abundance of phosphorus. More generally, inspection of the H I and S iv constraints shows 
that a similar conclusion applies for most metals. It is difficult to construct a multi U-zone 
model which will reproduce a better S iv to H I ratio without invoking a larger departure 
from solar metalicity for both higher and lower U models. 

For comparison, using the MF spectrum for the same analysis presented in Figure 7 
leads to the following conclusions: The absorber is characterized by 19.8 < log{NH) < 20.3 
and —2 < log(f/) < —1.5 and a minimum metalicity enhancement of 2-3 compared to solar. 
The metalicity enhancement is necessary since for the MF spectrum the strips allowed by 
the Siv and Hi constraints do not overlap, and most of the metals' lower limit curves lie 
somewhat above the allowed H l zone. 



4.1. The Phosphorus Over-abundance 

Previous claims for a large phosphorus over-abundance in PG 0946+301 are found in 
Junkkarinen et al. (1997). By comparing the derived column densities of Civ and P v and 
accounting for possible ionization differences, these authors derived (P/C)o946 ~ 6O(P/C)0. 
However, this assertion is crucially dependent on the assumption that the column densities 
deduced from converting the residual intensities to optical depth are a good approximation 
for the real ones. Even at the time, the possibility of saturation and covering factor was 
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addressed and Junkkarinen et al. acknowledged that: "If the Civ column density is under- 
estimated by a large factor and the P V column density is underestimated by a small factor, 
then the apparent over-abundance of phosphorus could be much less." Since there is now 
ample evidence for saturation and varying covering factor in the BALs of PG 0946-1-301, the 
phosphorus over-abundance claim based on the P v/C iv ratio cannot stand on its own. 

Our claim for a phosphorus over-abundance does not suffer from this weakness since 
we do not treat apparent iVjo„ (i.e., the template fitting values) as real A^jon- Instead we 
use Nion lower and upper limits, which are largely immune to saturation and covering factor 
effects. As can be seen in Figure 7 a comparison between the P V lower limit and the H I 
upper limit necessitates a large phosphorus over-abundance in order for the two constraints 
to overlap in parameter space. The exact value will depend somewhat on the shape of the 
ionizing continuum and on the spread of U values, but a phosphorus over-abundance by a 
factor of less than 10 relative to solar requires a rather contrived photoionization model. 
We also checked whether optically thick (bound-free) photoionization models could ease or 
eliminate the need for phosphorus over- abundance. Our conclusion is that these models do 
not improve on the optically thin results, where the details of this investigation are found in 
the Appendix. 

Finally, it must be noted that there remain major uncertainties in some of the atomic 
physics implemented in spectral simulation codes such as Cloudy. Among these are many 
of the dielectronic recombination coefficients, especially for 3rd and 4th row elements, such 
as P, S, and Ar (see Ferland et al. 1998). Until more reliable coefficients become known, 
any conclusions drawn regarding gas abundances of these elements should be done so with 
caution. 



5. SUMMARY 

We have presented an HST/STIS spectrum of PG 0946+301, which constitutes the 
highest quality UV data set of any BALQSO. These data are both high signal to noise and 
covers a very large number of BALs between 500 A and 1550 A in the rest frame. This 
provides for much tighter constraints on the ionization and abundances than previously 
possible. These data are available to the community for further analyses at: 
http://unix.cc.wmich.edu/korista/ftp/pg0946/spectra 

Our preliminary analysis of these data yield the following results: 

• Most of the 27 BAL troughs detected are saturated while not black (i.e., their shape 
is largely determined by partial covering effects). 

• The main outflow component is characterized by a total column density 20.3 < log(A^H) < 
20.8, ionization parameter — 1.7 < log(f/) < —1.2 and is consistent with solar abun- 
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dances. However, the metalicity can be up to ten times solar due to the large uncer- 
tainty in the H I BAL column density. These findings take into account the complica- 
tions posed by the saturation mentioned above. 

The abundance ratio of phosphorus to other metals is roughly ten times the solar ratio. 
Ionization models including a Hell ionization front cannot significantly decrease the 
necessary phosphorus over-abundance. 

In addition to the main outfiow component, a much higher ionization component ap- 
pears to be present. This wide range in ionization parameter is inferred from the de- 
tection of Ne VIII and supported by the constraints from other high ionization species. 

A new shallow high velocity trough appeared in these spectra that was not present in 
the 1992 HST/FOS spectra. This trough is evident only in the higher ionization lines, 
and the depth of this trough appears to increase with ionization potential of the ion. 
Significant changes in depth occurred in the main trough as well. 
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APPENDIX: CAN MODELS WITH IONIZATION FRONTS LESSEN THE 
REQUIRED PHOSPORUS OVER-ABUNDANCE? 

Due to the significance of a phosphorus over-abundance, we made efforts to check 
whether optically thick (bound-free) ionization models can explain the observed A^jon con- 
straints without invoking departure from solar metalicity ratios. The most promising set of 
such models are those that include a strong Hell ionization front (Voit, Weymann & Ko- 
rista 1993), where the incident spectrum is strongly suppressed beyond 54 eV, the ionization 
potential (IP) of Hell. Such an ionizing continuum may protect P v from photo-destruction 
(IP=65 eV), while still allowing for the photoionization of H I (IP=13.6 eV) and Si iv (IP=47 
eV) . In principle this can narrow the gap between the H I and P v curves. Another advantage 
of such a model is that it might explain the sharp dichotomy we observe between the low 
ionization species (Si IV, Nlll, Siv and Hi) and the higher ionization species (Civ, Nv, 
Oiv, Svi...), since the Hell ionization potential falls in between these groups. However, a 
possible problem is that Ar vil should also be protected by a He ll front (although to a lesser 
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extent because of its higher IP of 124 eV), but we do not detect it in the data, and due to 
the high oscillator strength of the Ar vil A585.75 singlet we obtained a strong constraint on 
the upper limit of its column density (< 10^^ cm~^). 

In order to test whether He II ionization front models can eliminate or lessen the need for 
phosphorus over-abundance, we constructed numorous CLOUDY models with such fronts. 
Here we give a brief summary of several representative models. We attempted to find solar- 
abundances CLOUDY models which simultaneously reproduce the inferred upper limits to 
the column densities for Hi, Siv and Arvil, and the lower limit for the Pv column density. 
Due to the low solar-abundance of phosphorus, CLOUDY models tend to produce higher 
ratios of H l/P v, S iv/P v and Ar vil/P v compared to the observed values, so that the goal 
of our modelling is to bring the upper limits of these more abundant species into agreement 
with the lower limits of P V. The observed column density ratio limits are quite constraining, 
since the ionization potentials of the other ions bracket the IP of P v (see above). 

We let CLOUDY find the optimal model (i.e., finding a minimum for the y^ between the 
observed and model-produced constraints, summed over the four constraints detailed above) 
by changing Nh and U . The main difficulty is due to the unknown shape of the incident 
ionizing continuum. We experimented with many continuum shapes and results from four 
representative models are shown in Table 3: 1) The canonical Mathews and Ferland (MF) 
AGN spectrum. 2) The F^j oc z/~^ ionizing spectrum described in § 4. 3) A very steep 
spectrum of F;^ oc v^^'^ above 1 Ry, which, although not realistic, comes closest to being 
consistent with the observed column density ratio limits 4) A more physically realistic version 
of the above, which takes into account the observed z/~^ spectral slope and connects to the 
MF spectrum at high and low energies. 

As can be seen in the table, the column density ratios of Ar vil and S iv to P v are 5-10 
times higher than observed. The main reason for this is that if a slab of gas contains a He ll 
front there are at least two ionization zones present: a zone with highly ionized gas in which 
He is completely ionized (the Hell Stromgren layer), and a second zone behind it in which 
He II becomes the dominant helium ion. The first highly ionized zone contains a significant 
amount of Arvil. The spectrum irradiating the second zone contains a Hell edge, and this 
zone contains virtually all of the P V. Due to the high U values of the models (necessary 
in order not to over produce Hi) the first zone is so thick that the ArVII column density 
constraint is already exceeded at a total hydrogen column density at which no detectable 
Pv is predicted. If the models try to lower U in order to produce less Arvil before the 
front, then too much Hi and SiV are produced. The end result is a compromise that keeps 
both Ar vil/P V and S iv/P v ratio about an order of magnitude above the observed values. 
Only a very contrived ionizing spectral shape can change this situation. We note that if the 
outflowing gas sees an ionizing spectrum which already includes a strong Hell front, it is 
possible to satisfy all four observational constraints. However, in that case we should see 
evidence for the "shield" that created the He ll front in the spectrum. Since such a "shield" 
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is very probably static, we should see deep absorption features corresponding to it in C IV, 
N V and O vi, close to the systemic velocity. Such features are not detected in the spectrum. 
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Table 3. photoionization models with He ii fronts 



# 


spectral shape 


log(iV,) 


^ log(f/) f 


Thi 


^ Th.u 


h Arvil 
Pv 


k Siv : 
Pv 


' iV,, - 


1^ 


MF 




21.46 


-1.1 


1.3 


190 


10 


4 


2 X 10^^ 


2"^ 


a = -2 




21.67 


-0.85 


0.9 


460 


10 


9 


1.4 X 10^^ 


3^ 


a = -3.5 




21.54 


0.00 


0.1 


330 


4 


4 


1.4 X 10^'^ 


4d 


a = -2, - 


-3.5 


21.52 


-0.76 


0.5 


290 


6 


6 


8 X 10^"^ 



^Mathews and Ferland (MF) AGN spectrum, with a 10 /im break. 
^F,y oc z/~^ for 1 Ry< hu <40 Ry and MF shape elsewhere. 
^F,y oc u^^-^ for 1 Ry< hu. 

'^Fy oc z/~^ for 1 Ry< hv < 1.5 Ry, F^ oc z/~-'^-^ for 1.5 Ry< hu < 30 Ry and MF shape 
elsewhere. 

''Optimal Nh for the slab in units of cm^^. 

^Optimal U value for the model. 

^Optical depth at the lyman edge. 

^Optical depth at the Hell lyman edge. 

'^Ratio of modeled A^Arvll/^P v normalized to the observed ratio. 

'Ratio of modeled A'siv/^P v normalized to the observed ratio. 

™Model value in units of cm~^. 
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